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basal  of  the  two.  Such  a  placement  of  the  more 
terrestrially  adapted  Ichthyostega  (76),  taken 
together  with  the  features  indicating  weight-bearing 
ability  (ventrally  facing  radial  and  ulnar  facets)  of 
the  very  primitive  ANSP  21350  (P),  would  sug¬ 
gest  a  scenario  of  rapid  early  terrestrialization 
rather  different  from  the  currently  predominant 
“aquatic  Devonian  tetrapods”  model. 
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A  Ferroelectric  Oxide  Made 
Directly  on  Silicon 
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Metal  oxide  semiconductor  field-effect  transistors,  formed  using  silicon  dioxide  and  silicon,  have 
undergone  four  decades  of  staggering  technological  advancement.  With  fundamental  limits  to  this 
technology  close  at  hand,  alternatives  to  silicon  dioxide  are  being  pursued  to  enable  new 
functionality  and  device  architectures.  We  achieved  ferroelectric  functionality  in  intimate  contact 
with  silicon  by  growing  coherently  strained  strontium  titanate  (SrTiOs)  films  via  oxide  molecular 
beam  epitaxy  in  direct  contact  with  silicon,  with  no  interfacial  silicon  dioxide.  We  observed 


ferroelectricity  in  these  ultrathin  SrTiOs  layers  by 
ferroelectric  nanodomains  created  in  SrTiOs  were 

For  decades,  semiconductor  device  design¬ 
ers  have  envisioned  numerous  devices 
using  ferroelectrics  in  combination  with 
semiconductors.  These  concepts  include  non¬ 
volatile  memories  (7,  2),  “smart”  transistors  that 
can  be  used  as  temperature  or  pressure  sensors 
(3),  and  ferroelectric  field-effect  transistors  whose 
logic  states  require  no  power  to  maintain  (4,  5). 
Missing,  however,  has  been  the  ability  to  in¬ 
tegrate  ferroelectrics  directly  with  mainstream 
semiconductors.  Our  work  bridges  this  gap, 
demonstrating  ferroelectric  functionality  in  a 
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means  of  piezoresponse  force  microscopy.  Stable 
observed  at  temperatures  as  high  as  400  kelvin. 

SrTiOa  thin  film  grown  directly,  without  any  in¬ 
termediate  layers  and  free  of  reaction,  on  the 
workhorse  of  semiconductor  technology,  silicon. 

Recent  work  has  explored  ways  to  use  epi¬ 
taxial  strain  to  induce  or  enhance  ferroelectricity 
in  thin  films  (6-8).  Ferroelectric  responses  that 
are  distinct  from  naturally  occurring  bulk  ferro¬ 
electrics  have  also  been  obtained  through  the 
growth  of  nano-engineered  superlattices  contain¬ 
ing  different  dielectric  and  ferroelectric  phases 
(9-11).  In  all  of  the  above-referenced  studies,  the 
substrate  and  film  are  isostructural.  For  commen¬ 
surate  SrTiO3/(001)  Si,  the  interface  is  structur¬ 
ally  far  more  complex,  connecting  a  diamond 
stmcture  (silicon)  with  a  perovskite  (SrTiOs)  (Fig. 
1).  The  high  reactivity  of  silicon  with  many  ele¬ 
ments  and  their  oxides  (72)  presents  a  formidable 
challenge  to  the  integration  of  functional  oxides 
with  silicon,  as  does  the  tendency  of  a  pristine 
silicon  surface  to  rapidly  form  its  own  oxide. 

Using  molecular  beam  epitaxy  (MBE),  we 
have  deposited  epitaxial  SrTiOs  films  on  (001)  Si 
substrates  via  a  kinetically  controlled  growth 
process  (73, 14)  (fig.  SI),  which  synchrotron  dif¬ 
fraction  measurements  reveal  to  be  commensu- 
rately  strained  up  to  a  thickness  of  ~24  A.  The 
growth  method  used  differs  substantially  from 
those  reported  previously  (15-1 7)  for  the  epitax¬ 


ial  integration  of  SrTi03  with  silicon  (73)  (figs. 
SI  and  S2).  Although  bulk  SrTi03  is  not  fer¬ 
roelectric  at  any  temperature,  the  large  compres¬ 
sive  strain  (~1.7%)  induced  on  commensurately 
strained  SrTi03/Si  is  predicted  to  result  in  fer¬ 
roelectricity  (73,  18)  with  a  Curie  temperature 
(Tc)  near  room  temperature  (fig.  S3  and  table  SI) 
and  an  out-of-plane  polarization  (6,  3, 19).  Films 
whose  thickness  exceeds  the  equilibrium  critical 
thickness  (20),  however,  are  unstable  to  relaxa¬ 
tion,  which  would  lower  the  transition  temper¬ 
ature  and  produce  nanoscale  heterogeneity. 

We  discuss  data  from  five  SrTi03  films  iden¬ 
tified  by  their  nominal  thickness  in  molecular 
layers  (ML):  5  ML,  6  ML,  8  ML,  10  ML,  and  20 
ML.  These  SrTi03  films  were  grown  on  (001)  Si 
substrates  by  MBL,  in  layers  of  one  to  a  few 
molecular  strata  at  a  time,  until  the  desired  thick¬ 
ness  was  reached  (73).  The  silicon  substrates  used 
in  this  study  were  doped  with  n-type  phosphorus 
(1  X  10^^  to  5  X  10^^  phosphorus  atoms/cm^) 
having  a  resistivity  of  1  to  4  ohm-cm.  Lach  layer 
that  was  grown  involved  a  controlled  sequence  of 


Fig.  1.  Structure  of  the  SrTiOs/Si  interface,  written 
and  imaged  on  a  6  ML  SrTiOs/Si  sample  by  PFM. 
With  the  45°  in-plane  rotational  offset  between  the 
unit  cells  (15),  the  epitaxial  orientation  relation¬ 
ship  is  (001)  SrTiOs  //  (001)  Si  and  [110]  SrTiOs  // 
[100]  Si. 
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steps  {12,  14)  that  kinetically  suppress  the 
oxidation  of  the  substrate  and  reduee  the  tend- 
eney  of  the  film  to  form  islands  (27)  (fig.  S4). 

X-ray  dififraetion  reveals  the  stmetural  quality 
and  strain  relaxation  that  oeeurs  in  the  SrTiOs 
films  as  thiekness  is  inereased.  Roeking  eurves  in 
CO  [where  the  angle  of  ineidenee  (0)  between  the 


x-ray  beam  and  the  sample  is  roeked  while 
leaving  the  deteetor  position  (20)  fixed]  of  the 
out-of-plane  SrTi03  002  refleetion  are  shown  in 
Fig.  2A.  Eaeh  eurve  displays  an  intense  and 
narrow  eentral  peak  due  to  eoherently  strained 
SrTiOs  on  top  of  a  broad  baekground  peak  (75). 
The  height  of  the  sharp  eentral  peak  in  relation  to 
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Fig.  2.  The  strain  state  of  the  SrTiOs  films  revealed  by  x-ray  diffraction.  (A)  Rocking  curves  in  co  of 
the  out-of-plane  SrTiOs  002  reflection.  (B)  Reciprocal  space  map  of  the  202  SrTiOs  peak  for  the 
6  AAL  sample.  Note  the  beating  along  the  I  direction.  (C)  Off-axis  scans  through  the  202  SrTiOs  peak. 
(D)  Reciprocal  space  map  of  the  202  SrTiOs  peak  for  the  8  AAL  sample.  In  (B)  and  (D),  the  intensity 
increases  from  white  to  red  to  blue  to  black. 


the  background  on  this  log  intensity  scale  gives 
an  indication  of  the  fraction  of  the  SrTi03  film 
that  is  coherently  strained.  The  coherently 
strained  fraction  of  the  SrTi03  films  decreases 
as  the  film  thickness  is  increased.  The  full  width 
at  half  maximum  (FWHM)  of  the  5  ML  sample, 
0.012°,  is  representative  of  the  sharpness  of  the 
coherent  peaks  (fig.  S5). 

We  used  x-ray  diflfaction  to  determine  the  in¬ 
plane  strain  of  the  SrTi03  films  (75, 22).  Because 
the  out-of-plane  lattice  constant  of  SrTi03  is  dis¬ 
tinct  from  that  of  silicon,  the  in-plane  lattice 
constant  of  SrTi03  can  be  obtained  by  measuring 
an  off-axis  SrTi03  reflection  where  there  is  no 
overlap  with  a  substrate  peak.  Figure  2C  shows 
scans  made  through  the  SrTi03  202  peak  for  the 
5  ML,  6  ML,  8  ML,  and  20  ML  samples.  The 
sharp  peak  observed  at  /z  =  ^  =  2.00  Si  reciprocal 
lattice  units  (r.l.u.)  is  due  to  the  commensurate 
portion  of  the  SrTi03  films  with  in-plane  lattice 
constant  =  ^si/  =  3.840  A.  As  the  film  thick¬ 
ness  increases,  the  relative  integrated  intensity  of 
the  sharp  peak  decreases  while  that  of  a  broad 
peak  dXh  =  k<  2.00  Si  r.l.u.  increases.  The  plot 
clearly  shows  the  transition  from  mostly  com¬ 
mensurate  SrTi03  to  mostly  relaxed  SrTi03  as 
the  film  thickness  is  increased.  Reciprocal  space 
maps  of  the  SrTi03  202  peak  for  the  6  ML  and 
8  ML  samples  are  shown  in  Fig.  2,  B  and  D, 
respectively.  The  6  ML  sample  (Fig.  2B)  has  its 
diffracted  intensity  mostly  centered  dX  h  =  k  = 
2.00  Si  r.l.u.,  whereas  for  the  8  ML  sample  (Fig. 
2D),  more  spectral  weight  is  observed  at  lower 
values  oi  h  =  k  (i.e.,  at  larger  in-plane  lattice 
constants)  because  of  relaxation  of  the  SrTi03. 
The  reciprocal  space  map  for  the  8  ML  sample 
also  shows  how  the  spectral  weight  tails  off  to 
higher  I  with  smaller  /z  =  ^  as  strain  relaxation 
sets  in.  From  f  scans  made  across  the  coherent 
peak  dit  h  =  k  =  2.00  Si  r.l.u.,  we  find  that  the 
coherent  peak  occurs  at  f  ~  2.71  Si  r.l.u. 

To  check  for  ferroelectricity  in  these  strained 
SrTiO3/(001)  Si  films,  we  used  piezoresponse 
force  microscopy  (PFM),  a  technique  that  has 
been  demonstrated  on  ferroelectric  films  as  thin 
as  28  A  (23-26).  With  strain  relaxation  occurring 
for  SrTiO3/(001)  Si  film  thickness  as  small  as 
8  ML  (~32  A),  measurement  of  the  piezoelectric 


Fig.  3.  PFM  images  (1  |xm  by  1  [im)  of  a  4  x  4  pattern  of  domains  written  on  the  5  ML  SrTiOa/Si  sample  at  different  temperatures.  (A)  T  =  298  K. 
(B)  r  =  314  K.  (C)  r  =  323  K. 
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Fig.  4.  Temperature  dependence  of 
the  out-of-plane  lattice  constant  of 
SrTiOs  strained  commensurately  to 
the  underlying  silicon  substrate.  Theo¬ 
retical  prediction  from  thermodynamic 
analysis  is  shown,  along  with  experi¬ 
mental  data  obtained  from  x-ray  dif¬ 
fraction  measurements  of  the  5  AAL 
sample.  The  error  bars  reflect  the 
maximum  error  expected  considering 
systematic  errors  and  variation  in 
sample  thickness.  Some  hysteresis 
between  data  taken  on  heating  and 
cooling  is  evident.  Also  indicated  is  the 
transition  temperature  observed  from 
PFAA  measurements  for  the  5  AAL 
sample  (T^s  aal  PFAA). 

0  50  100  150  200  250  300  350  400 

Temperature  {K) 


response  of  sueh  thin  layers  is  ehallenging.  The 
large  ^33  eoeffieients  predieted  for  strained 
SrTiO3/(001)  Si  (fig.  S6),  however,  make  it  a 
reasonable  signal  to  probe.  Loeal  eleetrie  fields 
were  applied  aeross  the  SrTi03  layer  by  means  of 
a  biased,  eondueting  atomie  foree  probe,  and  the 
resulting  piezoeleetrie  response  of  the  strained 
SrTi03  film  was  subsequently  imaged  using  the 
same  probe  {13).  At  room  temperature,  we  found 
that  domains  of  both  polarities  eould  be  patterned 
on  the  5  ML,  6  ML,  8  ML,  and  10  ML  samples, 
but  not  on  the  20  ML  sample.  Figure  1  shows  a 
PFM  image  written  on  the  6  ML  sample  at  room 
temperature.  In  all  of  the  samples  that  exhibited 
ferroeleetrieity  via  piezoresponse,  a  preferred 
downward  polarization  was  observed  (fig.  S7). 
This  agrees  with  reported  x-ray  fine  stmeture 
measurements  (27)  and  indieates  that  strained 
SrTi03  films  on  (001)  Si  are  prepoled  in  their  as- 
grown  state.  The  laek  of  observable  ferroeleetrie¬ 
ity  in  the  20  ML  sample  is  eonsistent  with  x-ray 
diffraetion  measurements  showing  that  the  20 
ML  sample  is  mainly  relaxed  (Fig.  2C)  and  helps 
rule  out  other  possible  meehanisms,  sueh  as 
ehanges  in  surfaee  ehemistry,  for  the  observed 
piezoeleetrie  response  in  other  samples.  A 
retention  study  of  the  written  domains  was  also 
earned  out  for  the  6  ML  sample  at  room  tem¬ 
perature  {13)  (figs.  S8  and  S9).  The  domain 
pattern  was  observed  to  be  stable  over  a  72-hour 
period,  at  whieh  point  the  pattern  was  erased  by 
rastering  the  atomie-foree  probe  with  a  eonstant 
voltage  over  the  patterned  area. 

PFM  measurements  performed  as  a  funetion 
of  temperature  revealed  a  rather  sharp  phase  tran¬ 
sition,  above  whieh  ferroeleetrie  domains  are  un¬ 
stable.  Figure  3  shows  a  series  of  three  PFM 
measurements  made  on  the  5  ML  sample  at  dif¬ 
ferent  temperatures.  Eaeh  image  was  aequired 
~30  min  after  writing  a  4  x  4  array  of  square 
domains.  Although  some  features  that  are  asso- 
eiated  with  imperfeetions  on  the  sample  surfaee 
also  showed  up  on  these  images,  at  T  =  298  K 
eaeh  of  the  16  domains  eould  be  observed  (Fig. 


3A).  The  temperature  was  inereased,  and  at  T  = 
314  K  only  9  of  the  16  domains  eould  be  seen 
(Fig.  3B).  The  existenee  of  a  single  domain  on 
the  third  row  from  the  top  mles  out  a  variety  of 
possible  measurement  artifaets,  sueh  as  a  “wan¬ 
dering”  eantilever  resonanee  frequeney.  At  T  = 
323  K  (Fig.  3C)  or  at  higher  temperatures,  no 
stable  domains  eould  be  observed.  These  PFM 
measurements  provide  a  lower  bound  on  the 
paraeleetrie-to-ferroeleetrie  transition  tempera¬ 
ture  (Tc):  Tc,5ml>314  K. 

Measurements  performed  on  the  6  ML 
sample,  however,  show  that  ferroeleetrie  domains 
written  on  it  are  stable  at  even  higher  temper¬ 
atures:  Tc,6  ml  >  410  K  (fig.  SIO).  Sueh  temper¬ 
atures  are  mueh  higher  than  that  predieted  by 
thermodynamie  analysis  (fig.  S3).  The  theoretieal 
ealeulation  assumes  an  infinitely  thiek  SrTi03 
slab  with  eomplete  polarization  eharge  sereening 
and  with  a  uniform  biaxial  eompressive  strain 
equivalent  to  that  obtained  by  growing  eommen- 
surately  strained  SrTi03  on  (001)  Si.  By  leaving 
out  surfaee  effeets  sueh  as  struetural  and 
eleetronie  diseontinuities  and  the  possibility  of 
ineomplete  sereening  of  the  polarization  eharge, 
the  thermodynamie  analysis  does  not  take  into 
aeeount  the  finite  film  thiekness,  whieh  presum¬ 
ably  would  lead  to  a  substantially  redueed  tran¬ 
sition  temperature,  as  has  been  shown  for  the 
related  ferroeleetries  PbTi03  {28)  and  BaTi03 
{29).  Thus,  the  observed  experimental  results  in- 
dieate  a  substantially  higher  transition  temper¬ 
ature  than  that  predieted  by  theory.  In  the  ease  of 
a  metal  in  eontaet  with  a  ferroeleetrie  {30,  31), 
polarization  sereening  at  the  interfaee  has  been 
shown  to  enhanee  the  ferroeleetrie  Tq.  Sereening 
of  the  polarization  eharge  as  well  as  struetural 
and  eleetronie  diseontinuities  at  this  heteroepi- 
taxial  SrTi03/Si  interfaee,  not  eonsidered  in  the 
present  thermodynamie  analysis,  eould  play  a 
role  in  understanding  the  quantitative  differenees 
between  experiment  and  theory. 

As  an  independent  eheek  of  the  ferroeleetrie 
phase  transition,  temperature-dependent  x-ray 
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diffraetion  measurements  of  the  out-of-plane 
lattiee  eonstant  were  performed  on  the  5  ML 
sample  (Fig.  4)  to  sense  the  stmetural  transition 
(7,  8,  32)  that  should  eoineide  with  Tq.  The  aver¬ 
age  out-of-plane  lattiee  eonstant  was  extraeted 
from  seans  made  of  the  SrTi03  002  peak.  With 
the  in-plane  lattiee  eonstant  elamped  to  the  sili- 
eon  substrate  and  ehanging  only  by  the  thermal 
expansion  of  silieon,  whieh  is  mueh  smaller  than 
that  of  SrTi03,  the  out-of-plane  lattiee  eonstant 
should  eontinuously  expand  with  temperature  in 
the  absenee  of  a  stmetural  transition.  The  mea¬ 
sured  out-of-plane  lattiee  eonstant  with  tempera¬ 
ture  for  the  5  ML  sample  shows  a  elear  deviation 
from  what  is  expeeted  for  thermal  expansion  with 
the  in-plane  lattiee  eonstant  eonstrained  to  that  of 
silieon.  This  deviation  eoineides  in  temperature 
with  the  transition  temperature  observed  by  PFM. 
The  “kink”  feature  observed  in  the  out-of-plane 
lattiee  eonstant  with  temperature  (7,  8,  32)  is 
qualitatively  eonsistent  with  thermodynamie 
analysis  of  eommensurate  SrTi03/Si  undergoing 
a  ferroeleetrie  transition  (Fig.  4),  although  the 
agreement  with  Tq  is  likely  to  be  eoineidental. 

A  ferroeleetrie  in  direet  eontaet  with  silieon 
invites  hybrid  ferroeleetrie-semieonduetor  dev- 
iees  {1-5).  Although  the  low  or  almost  nonexistent 
eonduetion  band  offset  predieted  {33)  and  mea¬ 
sured  {34)  between  SrTi03  and  silieon  eould  lead 
to  praetieal  diffieulties  in  implementing  sueh  fer¬ 
roeleetrie  deviees,  it  has  been  proposed  that  this 
problem  ean  be  overeome  by  earefully  eonstmet- 
ing  the  interfaee  between  SrTi03  and  silieon 
{35,  36). 
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Anomalous  Fractionations  of  Sulfur 
Isotopes  During  Thermochemical 
Sulfate  Reduction 

Yumiko  Watanabe/*  James  Farquhar,^  Hiroshi  Ohmoto^ 

Anomalously  fractionated  sulfur  isotopes  in  many  sedimentary  rocks  older  than  2.4  billion  years 
have  been  widely  believed  to  be  the  products  of  ultraviolet  photolysis  of  volcanic  sulfur  dioxide  in 
an  anoxic  atmosphere.  Our  laboratory  experiments  have  revealed  that  reduced-sulfur  species 
produced  by  reactions  between  powders  of  amino  acids  and  sulfate  at  150°  to  200°C  possess 
anomalously  fractionated  sulfur  isotopes:  =  +0.1  to  +2.1  per  mil  and  =  -1.1  to  +1.1  per 

mil.  These  results  suggest  that  reactions  between  organic  matter  in  sediments  and  sulfate-rich 
hydrothermal  solutions  may  have  produced  anomalous  sulfur  isotope  signatures  in  some 
sedimentary  rocks.  If  so,  the  sulfur  isotope  record  of  sedimentary  rocks  may  be  linked  to  the 
biological  and  thermal  evolution  of  Earth  in  ways  different  than  previously  thought. 


Large  anomalous  fractionations  of  sulftir 
isotopes  (1-5)  are  present  in  many  sedi¬ 
mentary  rocks  older  than  2.4  billion  years 
and  are  virtually  absent  in  younger  rocks  (6-10). 
It  has  been  argued  that  sulfur-bearing  minerals 
[such  as  pyrite  (FeS2)  and  barite  (BaS04)]  in 
sedimentary  rocks  older  than  2.4  billion  years 
formed  from  native  sulfur  (S^)  and/or  sulfate 
(S04^“:  produced  by  ultraviolet  (UV)  pho¬ 

tolysis  of  volcanic  sulfrir  dioxide  (SO2)  in  an  02- 
poor  atmosphere  (d,  8,  10),  and  thus  that  the 
record  of  anomalously  fractionated  S  isotopes  is 
evidence  for  the  transition  from  an  anoxic  to  oxic 
atmosphere  about  2.4  billion  years  ago  (7,  11). 
These  arguments  have  been  based  on  the  as¬ 
sumption  that  the  only  processes  producing 
anomalously  fractionated  S  isotopes  for  both 
A^^S  and  A^^S  (2,  6)  are  photochemical  reactions 
involving  gaseous  S-bearing  species  [such  as  hy¬ 
drogen  sulfide  (H2S)  and  SO2]  (6,  12).  Labora¬ 
tory  experiments  performed  with  UV  photolysis 
of  SO2  under  an  02-free  condition  produced 
and  S04^“  with  large  anomalous  fractionations  of 
S  isotopes  (13).  A  theoretical  study  also  suggests 
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a  maximum  partial  pressure  of  oxygen  (PO2)  of 
~10~^  atm  in  order  for  the  UV  photolysis  of  SO2 
to  produce  and  S04^“  (14).  Here,  we  present 


experiments  showing  that  reactions  between 
powders  of  amino  acids  and  S04^“  can  also 
produce  anomalous  fractionations  of  S  isotopes. 

FeS2,  the  most  abundant  sulfide  mineral  in 
sedimentary  rocks,  forms  from  a  variety  of  reac¬ 
tions  involving  H2S  and  Fe  in  sediments  and 
solutions  (15).  Both  bacterial  sulfate  reduction 
(BSR)  and  thermochemical  sulfate  reduction 
(TSR)  using  S04^“  and  organic  matter  in  waters 
and  sediments  are  important  in  the  production 
ofH2S 

2(C  -  C* *)  +  S04^“  +  2H^  ^  2CO2  +  2C*  +  H2S 

(1) 

where  C  and  C*  refer,  respectively,  to  reactive  and 
nonreactive  (residual)  carbon  atoms  in  organic 
compounds.  This  reaction  typically  results  in  posi¬ 
tive  correlations  between  the  FeS2  and  C*  contents 
of  sedimentary  rocks  (15).  The  C/C*  ratio  and 
reaction  rate  vary  depending  on  the  type  and 
maturation  degree  of  organic  compounds  (for  ex¬ 
ample,  carbohydrates,  hydrocarbons,  amino  acids, 
bitumen,  and  type  1, 11,  and  III  kerogens).  BSR  is 
carried  out  by  sulfate-reducing  bacteria  (SRB), 
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Fig.  1.  Cumulative  amounts  of  H2S  produced  from  SOs^  or  504^  reductions  by  amino  acids  (alanine  or 
glycine).  Gray  symbols  represent  the  experiments  using  alanine  (A)  and  solid,  open,  cross,  or  bar  symbols 
represent  those  using  glycine  (G).  The  three-digit  numbers  (150,  160,  170,  and  200)  represent 
experiment  temperatures,  and  the  hyphenated  numbers  represent  the  run  numbers. 
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Methods 


Film  growth 

SrTiOs  films  were  grown  using  elemental  strontium  and  titanium  sources  and  molecular 
oxygen  for  oxygenation.  Throughout  each  growth  the  quality  of  the  films  was  extensively 
monitored  in  situ  using  reflection  high-energy  electron  diffraction  (RHEED).  The  silicon 
substrates  (3”  diameter;  n-type  (phosphorous);  resistivity  =  1-4  Q-cm;  (001)  oriented 
within  +0.1°)  were  cleaned  using  a  commercial  UV  ozone  cleaner  for  ~20  minutes  before 
loading  into  the  molecular  beam  epitaxy  (MBE)  chamber.  The  chamber  pressure  was 
below  2x10'^  Torr  at  all  times  except  for  when  oxygen  was  introduced.  The  substrate 
temperature  (Tsub)  was  measured  by  an  optical  pyrometer  at  temperatures  above  550  °C 
and  by  a  thermocouple  (with  an  accuracy  of  +  50  °C)  at  lower  temperatures.  During  each 
deposition  step  discussed  below,  the  substrate  was  continuously  rotated.  The  native 
surface  oxide  of  the  silicon  substrate  was  thermally  removed  in  situ  prior  to  film  growth 
via  a  strontium-assisted  deoxidation  process  (SI)  by  depositing  ~2  molecular  layers  (ML) 
of  strontium  metal  after  heating  the  substrate  to  600  °C.  The  Si02  layer  desorbed,  with 
strontium  acting  as  a  catalyst,  when  the  substrate  was  heated  to  -800  °C.  The  oxide 
removal  was  observed  by  the  transformation  of  the  amorphous  RHEED  pattern  to  one 
showing  characteristic  features  of  a  crystalline  (001)  Si  surface  (see  Fig.  SI  A). 

With  this  deoxidation  scheme,  a  fractional  molecular  layer  of  strontium  remains 
bonded  to  the  silicon  surface  once  the  Si02  is  removed.  Its  presence  is  manifested  in 
RHEED  as  a  superposition  of  a  2x  and  a  3x  surface  reconstruction  along  the  [110] 
azimuth  of  (001)  Si  when  the  substrate  is  cooled  to  600  °C  (see  Fig.  SIB).  This 
superposition  of  surface  reconstructions  arises  from  a  strontium  coverage  between  1/6 
ML  (in  which  case  only  a  3x  reconstruction  would  be  observed)  and  Vi  ML  (in  which 
case  only  a  2x  reconstruction  would  be  observed)  bonded  at  the  silicon  surface  (S2,S3). 
More  strontium  was  deposited  at  600  °C  until  a  clear  2x  pattern  along  the  [110]  azimuth 
of  (001)  Si  was  observed  corresponding  to  a  total  of  Vi  ML  of  strontium  at  the  surface. 

Strontium  and  titanium  sources  were  carefully  adjusted  during  a  calibration  film 
grown  just  prior  to  the  actual  film  such  that  the  strontium  and  titanium  fluxes  were 
precisely  matched  at  -1x10^^  atoms/cm^-s.  This  precise  flux  matching  was  achieved  by 


first  using  a  quartz  crystal  microbalance  to  obtain  starting  values  for  strontium  and 
titanium  fluxes.  Using  these  initial  fluxes,  the  calibration  sample  was  grown  by 

_n 

codepositing  strontium  and  titanium  in  an  oxygen  background  pressure  >2x10“  Torr  and 
a  substrate  temperature  of  -700  °C.  During  this  ealibration  film,  using  RHEED  as  a 
guide,  the  temperature  of  the  strontium  cell  or  the  eurrent  through  the  Ti  Ball  {S4)  was 
adjusted  such  that  characteristic  surface  reconstructions  due  to  excess  strontium  or  excess 
titanium  (see  Fig.  S2)  did  not  appear  for  a  long  period  of  SrTiOs  growth.  This  ealibration 
process  yielded  a  relative  differenee  between  strontium  and  titanium  fluxes  of  less  than 
0.3%. 

The  oxygen  flow  was  metieulously  eontrolled  using  a  piezoeleetrieally  eontrolled 
leak  valve  to  give  a  moleeular  oxygen  flow  of  approximately  2x10^^  molecules/s  injected 
at  a  30°  angle  to  the  substrate  from  a  1  cm  diameter  tube  that  is  placed  -22  cm  from  the 
substrate.  If  this  flow  were  sustained,  the  chamber  pressure  would  eventually  reach  a 
steady  state  background  pressure  of  -1.5x10'^  Torr. 

With  the  substrate  temperature  at  -300  °C,  oxygen,  strontium  and  titanium  were 
eodeposited  to  form  2.5  ML  of  SrTiOa.  During  the  growth  (which  lasted  -2.5  min.)  the 
ehamber  baekground  pressure  typically  increased  into  the  upper  10'*  Torr  range  before 
the  oxygen  flow  was  stopped.  In  this  as-grown  state,  the  SrTiOs  film  is  epitaxial, 
although  the  crystal  quality  was  improved  by  heating  up  in  vaeuum  (less  than  2x10'^ 
Torr)  to  -580  °C  for  -10  min.  (see  Fig.  SIC  &  Fig.  SID).  To  grow  a  thicker  film,  the 
sample  was  cooled  down  to  -330  °C  and  a  further  1  or  2  ML  of  SrTiOs  was  added 
followed  by  another  vaeuum  anneal.  This  proeess  was  eontinued  until  the  desired 
thiekness  was  reaehed.  The  SrTiOa  films  grown  in  this  manner  showed  a  erystalline 
diffraetion  pattern  in  RHEED  at  each  step  of  the  sequenced  growth  (Fig.  SI).  X-ray 
diffraetion  measurements  earried  out  on  these  SrTiOs/Si  films  showed  erystalline  quality 
eomparable  to  single  erystal  SrTiOs  substrates  (Fig.  S5).  By  contrast,  the  growth  process 
discussed  in  ref.  S5  proceeds  through  recrystallization  of  an  amorphous  layer.  The 
amorphous  initiation  of  the  growth  in  this  latter  process  is  believed  to  impede  the  chances 
of  obtaining  a  commensurately  strained  SrTiOs  film  on  (001)  Si  (S6).  A  third  method  of 
growing  SrTiO3/(001)  Si  (ST),  which  takes  a  high-temperature  growth  approach,  was  also 
found  to  result  in  films  that  were  not  commensurate  with  silieon  showing  an  elastic 


anomaly  due  to  differential  thermal  expansion  between  SrTiOa  and  silicon  (S8). 
Moreover,  this  growth  scheme  resulted  in  an  extended  Si02  layer  at  the  SrTiOs/Si 
interface  (S9).  Although  optical  Raman  measurements  have  shown  evidence  for 
symmetry  breaking  in  films  grown  by  this  high-temperature  process  (S10,S11),  the  signal 
observed  was  attributed  to  dislocations  and  local  strain  in  the  vicinity  of  defects  in  the 
crystal  (SIO)  and  therefore,  cannot  be  attributed  to  a  consequence  of  ferroelectricity 
(Sll). 

Piezo-force  Microscopy 

Scanning  probe  microscopy  was  used  to  observe  ferroelectric  domains.  A  commercial 
atomic  force  microscope  (AFM)  (Asylum  MFP-3D)  was  employed  using  Pt-coated 
(OMCL-AC240TM-W2)  or  diamond-coated  (CDT-FMR-10)  silicon  cantilevers.  Piezo¬ 
force  microscope  (PFM)  images  were  acquired  using  the  following  protocol:  A  voltage 
was  applied  to  the  tip:  V{t)=Vtip+Vac  cos(M),  where  o)  is  an  angular  frequency  chosen  to 
be  close  to  the  resonant  frequency  of  the  cantilever-sample  system  while  in  contact 
mode.  The  ac  deflection  of  the  cantilever  was  measured  optically  and  detected  using  a 
lock-in  amplifier.  The  phase  of  the  lock-in  amplifier  was  calibrated  before  each 
measurement  such  that  a  positive  signal  of  the  in-phase  (X)  channel  was  observed  for 
sufficiently  large  Vtip  bias.  All  signals  acquired  were  taken  from  the  X  channel  in  this 
fashion.  To  read  ferroelectric  domains  Ftip  was  set  to  0  V. 

Ferroelectric  writing  was  achieved  using  the  following  protocol:  A  voltage 
Ftip(x,y)  was  specified  that  corresponds  to  the  image  one  intends  to  write  (e.g.,  atomic 
structure  of  SrTiOs/Si  interface  or  an  array  of  squares).  The  tip  was  held  fixed  while  the 
sample  position  was  scanned  in  a  raster  fashion  such  that  (x(0,y(0)  sweeps  out  the  entire 
area  to  be  written.  A  voltage  Vtip(.^(0>y(0)  was  applied  as  the  sample  was  rastered  to 
produce  the  desired  domain  structure. 


Thermodynamic  Analysis 


Fig.  S3  presents  a  phase  stability  diagram  for  (001)p  oriented  SrTiOa  thin  films  as  a 
function  of  the  in-plane  biaxial  strain  eo,  where  the  subscript  p  refers  to  the  pseudocubic 
Miller  index.  The  biaxial  strain  is  the  average  in-plane  strain  along  [100]p  and  [010]p 

axes  of  the  SrTiOa  film  from  its  underlying  substrate  and  eo>0  means  that  the  film  is 
under  tensile  strain. 

All  the  possible  stable  phases  in  the  biaxially  strained  (001);,  SrTiOs  films  shown 
in  Fig.  S3  are  listed  in  Table  SI.  In  the  table,  the  letters,  ‘T’  and  ‘O’,  indicate 
crystallographically  ‘tetragonal’  and  ‘orthorhombic’  symmetries,  respectively.  The 
superscripts,  ‘F’,  ‘S’  or  ‘F’,  distinguish  paraelectric,  antiferrodistortive  structural  (or 
ferroelastic),  and  ferroelectric  phases.  The  vectors,  p  =  (pi , F2 » Fa )  q  =  (q^,q2, ^3 ), 

are  the  spontaneous  polarization  and  the  linear  oxygen  displacement  corresponding  to 
simultaneous  out-of-phase  rotations  of  oxygen  octahedra  around  one  of  their  four-fold 
symmetry  axes.  They  are  the  order  parameters  to  describe  the  ferroelectric  transition  and 
the  antiferrodistortive  phase  transition.  Fig.  S3  is  obtained  using  thermodynamic  analysis. 
The  stable  phase  at  a  given  temperature  and  constraint  strain  minimizes  the  total  energy 
of  the  films.  Additional  details  maybe  found  in  ref.  S12. 

The  piezoelectric  coefficient  J33  (see  Fig.  S6)  was  also  calculated  using  its 
definition  of  J33  =  de^^jdE^ ,  where  is  the  applied  electric  field  component  on  the 
SrTi03  film  along  the  out-of-plane  direction,  ^33  is  the  strain  of  the  film  along  the  same 
direction  and  is  related  to  the  polarization  p  =  (0,0,  )  and  the  antiferrodistortive 

structural  order  parameter  q  =  (0,0,  q^ )  of  the  tetragonal  c-phase  (phase  Tj^  or 

shown  in  Fig.  S3)  by  ^33  =  pl{Q^^+2c^2Qnlcn)+ql{^n+^c^i^nlcn)-^c^2eJc^^  , 
where  Cij,  Qy  and  Ay  are  the  elastic  stiffness,  electrostrictive  coefficients  and  the  linear- 
quadratic  coupling  coefficients  between  the  strain  and  the  structural  order  parameter, 
respectively,  ps  and  q^  are  obtained  by  taking  the  minimum  of  the  total  energy  density  of 
f ^otal  f bulk  f elastic  ,  ^3 ,  Cq  )-  £'3P3  .  The  bulk  energy  density  and  the 

elastic  energy  density  of  the  film  are  given  in  ref.  SI 2. 


It  is  seen  that  the  piezoelectric  coefficient  ^33  (Fig.  S6)  increases  with  increasing 
temperature.  The  sharp  change  around  240  K  is  associated  with  the  disappearance  of  the 
AFD  structural  transition,  i.e.  ^3=0.  ^33  reaches  a  maximum  at  the  ferroelectric  transition 
temperature. 


X-ray  Diffraction  Measurements 

The  rocking  curves  through  the  SrTi03  002  reflection  presented  in  Fig.  2A  contain  both  a 
narrow  central  specular  component  and  a  broad  diffuse  scattering  component.  While  the 
full  width  at  half  maximum  of  the  specular  component  reflects  the  atomic  order 
coherency  length,  the  diffuse  intensity  describes  longer-range  correlated  structure  or  in 
some  cases  periodic  (coherent)  order.  For  the  20  ML  film,  the  rocking  curve  diffuse 
scattering  satellite  features  relate  to  long-range  lateral  domain  ordering  within  the  ordered 
oxide  film.  It  is  postulated  that  this  domain  ordering  forms  due  to  the  silicon  substrate 
surface  step  morphology  (S13)  as  the  SrTi03  fully  latches  to  the  underlying  Si  (001) 
surface. 

As  the  SrTi03  film  grows  45°  rotated  with  respect  to  the  in-plane  silicon  lattice 
(S14),  SrTiO3[100]//Si[l  10],  for  a  film  completely  clamped  to  the  silicon  substrate  the 
position  of  the  SrTi03  200  peak  occurs  at  the  exact  position  of  the  Si  220  diffraction, 
making  it  hard  to  distinguish  the  film  peak.  Since  the  1 10  reflection  is  forbidden  for 
silicon,  the  SrTi03 100  reflection  should  appear  background-free,  providing  a  means  to 
extract  the  in-plane  lattice  constant  of  the  SrTi03  film.  In  practice,  the  crystal  truncation 
rods  that  connect  the  allowed  111  and  111  silicon  reflections  pass  through  the  Si  1 10 
point  in  reciprocal  space  (or  close  to  the  Si  1 10  point  if  the  substrate  has  a  miscut), 
complicating  the  interpretation  of  the  data  at  the  SrTiOs  100  reflection.  Therefore,  to 
unambiguously  determine  the  in-plane  strain  of  the  SrTi03  films,  it  is  necessary  to 
measure  an  off-axis  SrTi03  reflection  (SI 5)  where  there  is  no  overlap  with  a  substrate 
peak  due  to  the  distinct  out-of-plane  lattice  constants  of  SrTiOs  and  silicon. 


Bistability  of  Ferroelectric  Polarization 


Fig.  S7A  shows  a  PFM  image  of  a  4x4  array  of  square  domains  written  on  a  pristine 
surface  (not  previously  scanned)  of  a  5  ML  thick  SrTiOs/Si  sample  using  a  tip  voltage  of 
Ftip  =  2  V  with  respect  to  the  silicon  substrate,  which  is  grounded.  A  similar  experiment 
performed  using  a  negative  bias  Ftip  =  -2  V  (Fig.  S7B)  shows  that  the  writing  of  stable 
domains  with  both  polarities  is  possible  on  the  same  area.  A  line  cut  across  both  images 
reveals  a  net  positive  piezoresponse  that  is  observed  in  the  areas  that  were  not  poled  in 
either  direction  (Fig.  S7C  &  Fig.  S7D).  This  positive  background  response  indicates  that 
the  as-grown  film  is  pre-poled  downward  (i.e.,  positive  charge  density  at  the  SrTiOs/Si 
interface).  This  result  is  consistent  with  ab  initio  calculations  as  well  as  x-ray  fine 
structure  measurements  (SI  6). 

Retention  Study 

A  pattern  consisting  of  four  squares  was  written  with  ytip=  ^  V  on  the  6  ML  thick 
SrTiOs/Si  sample  at  room  temperature  with  the  background  written  at  ytip=  +4  V.  The 
pattern  was  imaged  approximately  20  times  at  regular  intervals.  Fig.  S8  shows  some  of 
the  PFM  images  obtained  during  this  72-hour  period.  The  image  is  essentially  unchanged, 
except  for  a  small  lateral  shift  due  to  thermal  drift  of  the  AFM.  We  find  the  domain 
pattern  to  be  stable  over  a  72-hour  period.  There  is  a  reduction  in  contrast  due  to  changes 
in  the  resonant  frequency  of  the  cantilever;  the  spatial  resolution  is  also  apparently 
decreased,  again  due  to  wearing  of  the  AFM  tip  as  can  be  seen  from  the  topography 
images  (Fig.  S9)  obtained  simultaneously  with  the  PFM  images  of  Fig.  S8A.  Subsequent 
images  on  fresh  areas  showed  a  resolution  comparable  to  that  of  the  image  at  71  hr.  of 
Fig.  S8A. 

The  domain  pattern  was  subsequently  erased  by  applying  a  constant  Vtip=  ^  V  to 
the  same  area.  The  PFM  image  (Fig.  S8B),  taken  simultaneously  during  the  erase 
operation,  shows  a  slight  remnant  of  the  original  image.  Its  existence  arises  from  the  fact 
that  the  advancing  line  scan  has  only  erased  half  of  the  image  in  the  vicinity  of  the  probe. 
Subsequent  erase  scans  do  not  show  any  trace  of  the  original  domain  pattern,  nor  do  PFM 
images  taken  at  zero  bias  (ytip=  0  V,  Fig.  S8C). 


E 


Fig.  SI.  RHEED  images  at  different  stages  of  the  growth  of  SrTiOs/COOl)  Si.  (A)  After 
strontium-assisted  deoxidation  at  -800  °C  along  the  [110]  azimuth  of  (001)  Si  showing  a 
2x  reconstruction.  (B)  Superposition  of  2x  and  3x  surface  reconstructions  observed  as  the 
substrate  is  cooled  to  -600  °C.  Pink  and  orange  arrows  indicate  the  streaks  due  to  the  2x 
and  3x  reconstructions,  respectively.  (C)  At  rsub~300  °C,  the  as-grown  2.5  ML  thick 
SrTiOs/Si  film.  (D)  The  2.5  ML  thick  SrTiOs/Si  film  surface  after  heating  to  Tsub-SSO  °C. 
(E)  The  5  ML  thick  SrTiOa/Si  film  at  7’sub~300  °C.  (F)  Einal  surface  of  the  5  ML  thick 
SrTiOs/Si  film  after  heating  to  Tsub-SSO  °C. 


Fig.  S2.  RHEED  patterns  observed  from  a  stoichiometric  SrTiOs  surface  (top)  and 
characteristic  surface  reconstructions  observed  that  point  to  excess  strontium  or  titanium 
during  SrTiOs  growth  (bottom).  (A)  Along  the  [110]  azimuth  of  SrTiOs.  Orange  arrows 
indicate  extra  streaks  that  lead  to  a  2x  reconstruction  due  to  excess  strontium.  (B)  Along 
the  [100]  azimuth  and  (C)  along  the  [210]  azimuth  of  SrTiOs.  Blue  arrows  indicate  extra 
streaks  that  lead  to  2x  reconstructions  due  to  excess  titanium. 


SrTi03/Si 


In-plane  biaxial  strain 

Fig.  S3.  Phase  diagram  for  SrTiOa  under  varying  amounts  of  biaxial  strain. 
Commensurate  growth  on  silicon  corresponds  to  -1.7%  compressive  strain  at  room 
temperature  (indicated  at  left).  The  predicted  phase  transition  is  predicted  to  occur  just 
below  room  temperature. 


Fig.  S4.  Cross-sectional  annular  dark  field  (ADF)  scanning  transmission  electron 
microscopy  (STEM)  images  of  SrTiOs  films  grown  on  (001)  Si.  (A)  The  image  of  the 
nominally  5  ML  thick  SrTiOa  film  shows  an  average  film  thickness  of  6.5  unit  cells 
suggesting  a  non  uni  form  SrTiOs  coverage.  In  the  area  marked  with  an  arrow  more  than 
one  SrTiOa  grain  is  imaged  in  projection.  (B)  Lower  and  (C)  higher  magnification  images 
of  the  nominally  20  ML  thick  SrTiOs  film  confirming  the  absence  of  an  extended 
interface  layer.  The  ADF  images  were  recorded  on  a  200  kV  FEI  Tecnai  F20-ST  STEM 

O 

with  a  minimum  probe  size  of  1.6  A,  a  convergence  semiangle  of  (9.5+0.5)  mrad  and  an 
inner  detector  angle  of  65  mrad.  To  increase  signal  to  noise  and  average  out  the  scan 
noise  for  the  higher  magnification  images  (A)  and  (C),  10  and  6  successive  images,  each 
recorded  at  10  and  8  microseconds  per  pixel,  were  cross-correlated  and  averaged. 
Subsequently,  the  1024x1024  pixel  image  was  rebinned  to  512x512  pixels. 
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Fig.  S5.  X-ray  diffraction  rocking  curves  in  (o  of  the  out-of-plane  SrTiOa  002  reflection 
of  single  crystal  SrTiOa  substrates  and  of  the  10  ML  thick  SrTiOa  film  grown  on  (001)  Si. 
The  narrow  rocking  curve  of  the  SrTiOs/Si  film  indicates  its  excellent  crystalline  quality, 
eomparable  to  single  erystal  SrTiOa  substrates. 
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Fig.  S6.  Predicted  piezoelectric  coefficient  as  a  function  of  temperature  for 
commensurate  SrTiOs/COOl)  Si.  See  section  on  Thermodynamic  Analysis  for  details  of  the 
calculation  leading  to  this  prediction. 
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Fig.  S7.  Bistable  piezoelectric  response.  (A)  PFM  image  of  a  4x4  array  of  square  positive 
domains,  written  with  Vtip  =  +2  V  and  imaged  at  Vnp  =  0  V.  (B)  Same  as  (A)  except  that 
the  writing  voltage  Vtip  =  -2  V.  (C)  Linecut,  shown  above  as  a  dashed  line,  indicating  the 
profile  of  the  piezoresponse  for  the  case  where  Vtip  =  +2  V.  (D)  Same  as  (C)  except  Vup  = 
-2  V.  The  positive  background  piezoresponse,  in  the  absence  of  poling,  indicates  a 
preferred  downward  orientation  of  the  polarization  of  the  as-grown  film. 


Fig.  S8.  Retention  of  ferroelectric  domains.  (A)  Images  taken  at  different  times  during 
the  72-hour  period  of  a  2  pm  x  2  pm  area  of  the  6  ML  thick  SrTiOs/Si  sample  that  was 
patterned  with  four  square  domains.  (B)  PFM  image  obtained  as  the  domains  are  being 
erased  with  Vtip  =  ^  V.  Note  the  change  in  the  color  scale.  (C)  PFM  image  taken  after 
erasure  showing  no  trace  of  the  original  domain  pattern. 


Fig.  S9.  2  [am  X  2  |am  AFM  topography  images  of  the  6  ML  thick  SrTiOs/Si  sample 
taken  simultaneously  with  the  PFM  images  of  Fig.  S8A  during  a  period  of  72  hours. 
Images  show  a  decrease  in  spatial  resolution  with  time  due  to  wearing  of  the  tip. 
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Fig.  SIO.  1  |j.m  X  1  |j,m  PFM  image  of  four  square  domain  patterns  written  and  imaged  on 
the  6  ML  thick  SrTiOa/Si  sample  at  7=300  K  and  at  7=410  K. 


Table  SI.  Phase  definitions  in  biaxially- strained  (001)^  SrTiOa  films. 


Designation 

Point 

Group 

Symmetry 

Order  Parameter 

(p;q) 

(with  respect  to  the 
pseudocubic  cell) 

Equivalent  Domain  Variants 

(with  respect  to  the  pseudocubic  cell) 

jp 

4/mmm 

(0,0,0;  0,0,0)p 

4mm 

(0,0,0;  0,0,^3);, 

(0,0,0;0,0,-^3);, 

Of 

mm2 

(0,0,0;  ^1,0,0)^ 

(0,0,0;  -^1,0,0);,,  (0,0,0;  0,^i,0);„ 

(0,0,0;  0,-qu0)p 

of 

mm2 

(0,0,0,  ^1,^1, 0)jp 

(0,0,0,  -q\,q\,0)p,  (0,0,0,  q\,q\,0)p, 

(0,0,0;  -qi,-qi,0)p 

rpF 

4mm 

(0,o,p3;  0,0,0);, 

(0,0,-p3;  0,0,0];, 

Of 

mm2 

(pi,0,0;  0,0,0)p 

(-pi,0,0;  0,0,0);,,  (0,pi,0;  0,0,0)p, 

(0,-pi,0;  0,0,0)p 

^2 

mm2 

(pi,pi,0;  0,0,0)p 

(-pi,pi,0;  0,0,0)p,  (pi,-pi,0;  0,0,0)p, 

(-pi,-pi,0;  0,0,0)p 

rpF 

■*■2 

4mm 

(0,0, Pi;  0,0,q3)p 

(0,0,-p3;  0,0,q3)p,  (0,0,P3;  0,0,-q3)p, 

i0,0,-p3;  0,0,-q3)p 

Of 

mm2 

(pi, 0,0;  0,0,q3)p 

(-pi,0,0;  0,0,q3)p,  (pi,0,0;  0,0,-q3)p, 

(-pi,0,0;  0,0,-q3)p,  (0,pi,0;  0,0,q3)p, 

(0,-pu0;  0,0,q3)p,  (0,pi,0;  0,0,-q3)p, 

(0,-pi,0;  0,0,-q3)p 

Of 

mm2 

(puPuO;  0,0,q3)p 

(puPuO;  0,0,-q3)p,  (-puPuO;  0,0,q3)p, 

(-puPuO;  0,0,-q3)p,  (pu-puO;  0,0,q3)p, 

(pu-p\,0;  0,0,-q3)p,  (-pu-p\,0;  0,0,q3)p, 

(-pu-p\,0;  0,0,-q3)p 

Of 

mm2 

(puPuO;  ququ0)p 

(-pi,-pi,0;  q\,q\,0)p,  {p\,p\,0;  -qi,-qi,0)p, 

(-pi,-pi,0;  -qi,-qi,0)p,  (pi,-pi,0;  qi,-qi,0)p, 

(pu-p\,0;  -qu-q\,0)p,  (-pi,pi,0;  -quq\,0)p, 

(-puPuO;  qu-q\,0)p 
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